Stability of trapped fermionic gases with attractive interactions 
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We present a unified overview, from the mean-field to the unitarity regime, of the stability of a 
trapped Fermi gas with short range attractive interactions. Unlike in a system of bosons, a Fermi 
gas is always stable in these regimes, no matter how large the particle number. However, when the 
interparticle spacing becomes comparable to the range of the interatomic interactions, instability is 
not precluded. 



I. INTRODUCTION 

The freedom to tune the interparticle interactions in 
dilute Fermi and Bose gases using Feshbach resonances 
raises the question of the stability of trapped atomic gases 
with attractive interactions against collapse. Weakly 
interacting Bose gases with attractive interactions be- 
come mechanically unstable once the particle number 
exceeds a critical value, of order the ratio of the os- 
cillator length to the magnitude of the s-wave scat- 
tering length 

B II II S- The question is whether 
Fermi systems are always stabilized, and by what mech- 
anisms, or can they ever be similarly unstable? Our 
aim in this paper is to present an overview of the prob- 
lem in the different parameter regimes, putting together 
and extending the various arguments that have been 
iven for the stabil ity (an d instabilit y) o f Fermi systems 

i H H 113, m mEl El El El llTl within a unified 

framework. 

At first, one might think that the Fermi energy of the 
atoms would be sufficient to support a trapped Fermi gas 
against collapse, in the way, for example, that electron 
degeneracy pressure supports white dwarfs against grav- 
itational collapse. In a white dwarf the gravitational en- 
ergy scales as the inverse of the system radius, 7?., while 
the Fermi energy scales as 7?,"^ non-relativistically; in- 
stability sets in only in the relativistic limit, when the 
Fermi energy begins to scale as TZ~^. In a trapped 
Fermi gas with short range s-wave interactions between 
different components, the Fermi energy similarly scales 
as TZ~^, while the interaction energy scales, within mean- 
field theory, as TZ^^. Thus the Fermi energy is not, in 
itself, adequate to prevent an instability. 

In the Fermi problem with low energy s-wave interac- 
tions, one must distinguish several regimes. The first is 
that of the weak interactions where the s-wave scattering 
length, a, is small compared with the interparticle spac- 
ing, and mean-field theory is valid. Calculating within 
mean-field theory one finds an instability at high density 
^, 15, 16, 17J, kp\a\ > 7r/2, or large particle number 
N^^^^0.6d/\a\ |a,l3> where kp is the Fermi momentum, 
and d = ^/h/niLo is the oscillator length, with lo the os- 
cillator frequency (assumed isotropic); however, at such 
densities the interparticle spacing is comparable to the 
scattering length, and the mean-field calculation becomes 
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FIG. 1: Regions of stability in the kpro, kF\a\ plane of a two 
component Fermi gas with short range interactions. Here kp 
is the Fermi momentum, a the s-wave scattering length, and 
ro the range of the interactomic potential. 



invalid. Indeed, in the mean-field regime (fecial << 1) 
the system is always stable. The second - "unitary" - 
regime is when the scattering length is large compared 
with the interparticle spacing, itself large compared with 
the range of the interparticle potential. Here the sys- 
tem is arain always stable, since the total energy remains 
positive H milaESElElElEl- The third regime 
is when the interparticle spacing becomes comparable to 
the range, tq, of the interparticle potential. In this regime 
the system can in principle become unstable. Figure ^ 
schematically illustrates these regimes in the kpro, kp\a\ 
plane. 

We first study the stability of the Fermi gas in mean 
field theory, working in the local-density approximation 
in which the gas is assumed at each point to be locally 
homogeneous [lq| , and derive analytically the mean-field 
limits on the density and particle number. As we see, 
the two instability conditions, kp\a\ > 7r/2, and number 
N^/^>OM/\a\, are equivalent. We next discuss the uni- 
tary regime. By means of a simple interpolation formula 
for the density dependence of the energy, we show that 
the gas should be stable, with no limit in the number of 
particles in the system. We show, by using the Hartree- 
Fock approximation as a rigorous upper bound to the 
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total energy of a system with an attractive square well 
interparticle interaction, how stability can break down in 
the high density regime. 



II. STABILITY IN MEAN-FIELD 

We consider a Fermi gas with equal numbers of two 
components (denoted by | and I) in equilibrium at 
zero temperature in an isotropic harmonic trap, V{r) = 
^muP'r^ . The energy of the system in the local density 
approximation is 



E = 



10m 



(37r2)2/3 / dv n^/ 



-—mu! 
2 



dr n + g / dr n^ni 



(1) 



where g = AttTt?' a / m is the coupling constant, and n = 
2n| — 2ni is the total density of particles. The short 
range s-wave interactions between particles of the same 
species essentially vanish due to the antisymmetry of the 
many-fermion wave function. 

The equilibrium density minimizes the total energy 
with respect to the density, at fixed particle number, 
N = J dr n{r). Taking the constraint of fixed total num- 
ber of particles into account via a chemical potential, 
fi = 11^ ^ fj-l, we find that the first functional derivative 
of the energy is 

SE -fiSN^ [dr( (37r2)2/3^n2/3 + Irmj^^ 
J \ 2m 2 

+ ln-^^^^n^O, (2) 
which yields the density-position relation, 

(3^2)2/3^^2/3 _^i^^2^2 5^^ (3) 

^ ' 2m 2 2 ^ ^ ' 

It is most convenient to regard this equation as deter- 
mining in terms of n: 

^- ((3^2^)2/3 - 4^|a|n) + R^{a, N), (4) 



where d = y^Ti/muj is the oscillator length, and the radius 
of the cloud is i? = y^2/I/mc<j2. Equation Q can be 
inverted either graphically, or by simply solving the cubic 
equation explicitly. The results are shown in Fig. |5J). 
Note that when a < 0, r(n)^ has a minimum as a function 
of density at 
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FIG. 2: Distance from center of trap, squared, as a function 
of the density, n, Eq. Q i for the value R — d. For large nega- 
tive a the curve does not intersect the n-axis, and the system 
becomes unstable. The dashed curve shows the characteris- 
tic radius squared vs. density calculated with the effective 
scattering length, Eq. (IIH . 



In order for the density profile to be stable it must 
minimize, not maximize, the energy functional, i.e., the 
second variation of £' — fiN: 
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dr (37r^) 
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must be positive for all Sn. This condition requires. 



kFir)\a\ < 



1.57, 



(7) 



(8) 



at all positions, or equivalently n < 7r/24|ap (cf. 
Eq. ©). As we see, Eq. JHJ is the condition that the 
density profile corresponds to a minimum of the energy 
functional Note that we have derived this result in 
terms of the Thomas- Fermi profile including interactions. 
This argument is similar to that of The result © 
has also been obtained heuristically by comparing the re- 
pulsive force on an atom arising from the kinetic energy 
term with the attractive force due to the mean field in- 
teraction energy 0, and derived via a field-theoretical 
approach [TH- 

When the minimum of vs. n is at negative r^, the 
system is stable because then n is less than no, so that 
kp\a\ is always < 7r/2. On the other hand, when the min- 
imum is at positive values, the curve does not intersect 
the horizontal axis, and therefore there is no solution for 
the density at the center of the trap. (An unphysical so- 
lution for the density would have a hole of radius tq about 
the center of the trap.) However, in this case the den- 
sity equals no at r = tq, and vanishes as r —^ R. Then, 
A:F(0)|a| = 7r/2, and the solution is marginally stable, 
not actually minimizing the energy functional We 
conclude that the system is stable when rg < 0, and the 
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critical point is tq = 0, corresponding to 

12|a|2- 



(9) 



and then the density of the system at the center of the 
trap is given by Eq. JB)). 

When Eq. is satisfied the system reaches its max- 
imum number of particles in equilibrium, given by the 
area under the n vs r curve, Nmax — 47r n[r)r'^dr. 
Integrating with respect to the n, we find 



N„ 
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(10) 



where 
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1/2 

(1 — x)x^dx 
! 0.00823. 



As long as the local density approximation is valid, the 
condition (|1U|) is equivalent to the condition ||SJ|. How- 
ever, the use of mean field is no longer valid for fci?|a|;il. 
From IjlOl) . the condition for stable equilibrium is N^^^ < 
0.612d/|a|, which agrees with that found numerically in 
Ref. 0. 

Table ^ lists the number of particles in representa- 
tive experiments on ^Li gases with two-spin components. 
Since these experiments were not testing the maximum 
number of particles that could be trapped, they provide 
only lower bounds on the maximum number of particles. 
For small kp\a\ < 1, the experimental particle numbers 
are well within the bound H10|) . However fermionic sys- 
tems at larger fc_F|a|, exceeding tt/2, do remain stable. 
See Fig. |3| The argument leading to Eq. (|10|l breaks 
down at large fcF|a|, for then correlations beyond mean 
field become important. In addition, when n(0) — > uq, 
the local density approximation fails. A necessary con- 
dition for the local density approximation to be valid is 
that the density varies slowly over an interparticle spac- 
ing, or \d\nn/dr\ <C kp- From Eq. 10}, we thus find the 
requirement, for n(0) ~ no, that 1 — 2kF\a\/Tr ^ (|a|/d)2 
for the local density approximation to be valid, where we 
have taken to be bounded by the critical in Eq. J^. 
This condition clearly breaks down as n{0) hq, where 
dn/dr diverges, as one sees in Fig. |2 



III. NEAR THE UNITARY REGIME 

When kp\a\ becomes large perturbation methods fail, 
and the question of stability is more involved. In the 
unitarity regime (/cF|a| — > oo), the energy of the ground 
state is bounded and can be written in the scale- free form 
01111113 -Eo = (1 + 0)Efg, where /3 is universal con- 
stant —0.5. To approach the problem, we interpolate 
the energy between small \a\ and — oo by calculating in 
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FIG. 3: The maximum particle number vs. —a/d, as pre- 
dicted by Eq. IIUH . The nearly horizontal lines are con- 
tours of constant chemical potential in the stable region; here 
R/d = yj2^/hLu. Also shown are the particle numbers mea- 
sured in the five experiments in Table 1. Data points are 
labelled by the corresponding reference numbers. 



TABLE I: Experimental measurements of particle numbers in 
''Li gases. In the regime kF\a\ < 1, Eg. 1)10^ is satisfied. 



Expt. 


T/Tf 


-a/ao 




max 




m 


< 0.2 


3500 


9 X 10^ 


-10^2 


< 1 


[20] 


< 0.1 


~ 500 


6 X 10^ 


~ 10^ 


<< 1 


M 


~ 0.75 


15600 


1.3 - 9 X 10** 


~ 7 


>> 1 


m 


~ 0.15 


10^ 


8.0 X lO"* 


~ 1 


^ 2.5 


m 


~ 0.1 


10^ 


7.5 X W 


~ 


^ 7.4 



mean field theory with a density dependent effective scat- 
tering length, 



fteff 



1 — jkpa ' 



(11) 



for kpa < 0, and replacing the couphng g in Eq. ^ by 
the local coupling g^s = 47r?l^aefi/m. The choice 7 = 
— 20/97r/3 reproduces the known energies in the weakly 
interacting limit and in the unitary regime |2lj| . 

Repeating the above Thomas-Fermi calculation with 
this density-dependent scattering length, we find, first, 
the position as a function of density. 



k p 



47r|acff|n - '^kpla^ff 



d^ + R^. (12) 



The density profile is similar to that of the stable config- 
urations shown in Fig. 2. We show in that figure, as a 
dashed line, the density profile calculated with Eq. H12|) 
a = -0.5, and P = -0.5. 
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FIG. 4: The second variation oi E — fiN, Eq. (I13II . as a 
function of —l/kF\a\, calculated with the effective scattering 
length, Bcfi, Eq. llTTt . 



The second functional derivative of the energy with 
respect to n is, 



S^E - ^i6^N 



10 



(tt + 2kFaeS + -z-likpa, 



off') 



+ ^l\kpa,sf). (13) 

In the hmit 7A:i?|a| << 1, a^s ~ a, and Eq. 
predicts that the gas is stable for all /cF|a|, in dis- 
tinction to the prediction of Eq. ||SJ|. In the opposite 
limit, 7A:F|a| >> 1, locffl ^ %iTj3 /20kp, a constant, 
and {5^E - ii5^N) /{5n^) ^ (7r'^h^)/{mkF){l + (3/2). 
Hence the gas is stable for /? > —2, as can be seen in 



fiS"^ N) / {Sn)"^ as a 
6Li 



Fig. H a plot of imkF/'n-ri^)i5^E 
function of —l/kp\a\ for various /?. Experiments on 
find [HI fill IT^ . P ^ —0.5; for these values, the second vari- 
ation is always positive, predicting that the gas should 
be stable for all fc^laj, in agreement with experiment. 



IV. EVEN HIGHER DENSITY 

In the scale-free unitary regime the interparticle spac- 
ing is large compared with the range of the potential, 
and the system is stable. However, at higher densities 
when the interparticle spacing becomes comparable to 
the range of the potential, the system can in fact become 
unstable, as we see from the following model calculation. 
Let us assume that the interparticle interaction is an at- 
tractive square well potential with range ro and depth 
V{< 0). The energy calculated assuming a Hartree-Fock 
wave function is a rigorous upper bound to the exact 
total energy, E; thus 



E , Sh^kl 



< 



(14) 



TV - lOm ■ 3 
As n increases, the energy becomes unbounded below; 

the kinetic energy does not save the system from col- 
lapse. We note that a,t V — — ?i^7r^/4mrQ the two 
body-scattering amplitude has a resonance, the scatter- 
ing length diverges, and the interaction energy per parti- 
cle becomes Eint/N = —{5n/54)kproEFG- The instabil- 
ity in this model should not occur, however, for realistic 
interatomic potentials which have a repulsive core. 

We thank Randy Hulet and Zhenhua Yu for useful 
discussions. This work was supported in part by NSF 
Grants PHY03-55014 and PHY05-00914. 
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